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I.  Introduction 


The  bearingless  main  rotor  is  the  design  configuration  for  the  Army’s  next  generation  combat  heli¬ 
copter;  it  offers  design  simplicity  (fewer  parts),  weight  reduction,  better  maintenance,  and  more  control 
power  and  maneuverability,  A  bearingless  rotor  is  a  special  case  of  a  hingeless  rotor,  in  which  the  pitch 
bearing  as  well  as  the  flap  and  lag  hinges  are  eliminated  (see  Figure  1).  Pitch  control  from  the  pitch  link 
to  the  main  blade  is  transmitted  via  a  torsionally  stiff  torque  tube.  This  in  turn  twists  a  torsionally  soft 
flexbeam,  which  functions  effectively  as  a  pitch  bearing.  The  flexbeam  is  actually  a  major  component  of 
the  bearingless  rotor  design;  it  carries  the  centrifugal  load  and  allows  for  the  blade  flap,  lead  lag,  and  twist 
motions.  There  is,  however,  a  lack  of  vulnerability  information  for  the  flight-critical  flexbeam  element  at¬ 
tributable  to  ballistic  damage  mechanisms.  In  addition,  the  effects  of  flexbeam  ballistic  damage  on  the  rotor 
and  helicopter  system’s  performance  are  not  adequately  understood. 


The  objective  of  the  present  research  is  to  analytically  investigate  rotor  and  aircraft  performance  ef¬ 
fects  caused  by  flexbeam  damage.  The  present  investigation  was  performed  using  the  University  of  Maryland 
advanced  rotorcraft  code  (UMARC)  [1].  The  bearingless  main  rotor  (BMR)  system,  including  flexbeams, 
torque  tubes,  and  main  rotor  blades,  was  modeled  as  a  number  of  elastic  beam  finite  elements,  wherein  each 
beam  element  undergoes  flap  bending,  lag  bending,  elastic  twist,  and  axial  deflections.  Flexbeam  ballistic 
damage  was  simulated  as  changes  in  the  span- wise  distribution  of  the  mass,  bending  and  torsional  stiffness 
of  flexbeam  element. 


Results  are  first  calculated  for  an  advanced  five-bladed  BMR  with  a  baseline  (undamaged)  configura¬ 
tion.  Results  then  are  calculated  for  this  rotor  system  with  representative  levels  of  ballistic  damage  to  the 
hub  flexbeam.  The  effects  of  this  damage  on  this  helicopter’s  performance  are  determined  in  terms  of  blade 
modal  shapes  and  frequencies,  rotor  system  dynamic  response  and  loads. 


In  the  context  of  the  U.S.  Army  Research  Laboratory’s  (ARL)  process  structure  for  analyzing  combat 
system  vulnerability,  this  study  and  its  associated  engineering-based  methods  address  the  mapping  from 
Level  2,  the  Target  Component  Damage  State,  to  Level  3,  the  Target  Capability  State  (i.e.,  02,3  mapping); 
other  levels  include  1,  Threat/Target  Initial  Conditions  and  4,  Target  Combat  Utility.  See  Reference  2 
for  more  details.  A  significant  feature  of  this  process  is  that  at  each  level  (or  space),  distinct,  measurable 
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information  is  available,  defining  the  threat-target  encounter  and  vulnerability /lethality  outcome.  Here,  for 
example,  physical  and  structural  factors  defining  hub  flexbeam  damage  (Level  2)  are  mapped  via  engineering 
methods  into  parameters  that  define  the  rotor  and  helicopter  system’s  functional  capability  (Level  3);  all 
the  defining  terms  are  explicit  and  measurable  through  experimentation.  Application  of  these  and  other 
engineering  analysis  tools  to  the  vulnerability /lethality  process  structure  will  occur  largely  through  imple¬ 
mentation  in  the  degraded  states  vulnerability  methodology  for  level  02,3  mapping  now  being  developed  for 
aircraft  targets  at  ARL. 


II.  Nomenclature 


Cu, 


Ely 

Eh 

GJ 

h 

m 

rrio 

q 

R 

T 

u 

V 

V 
w 
W 
a 
C^s 
00 

0i8)01C 


Helicopter  weight  coefficient,  W/7rR?p{QR)^ . 

Blade  mean  chord,  ft. 

Flap-wise  bending  stiffness,  lb  —  in^. 

Chord- wise  bending  stiffness,  lb  —  m^. 

Effective  sectional  torsional  stiffness,  lb  —  in^. 

Length  of  the  zth  beam  element,  ft. 

Blade  section  mass,  slug. 

Reference  blade  section  mass,  slug. 

Number  of  blades. 

Global  displacement  vector. 

Blade  radius,  ft. 

Rotor  thrust,  lb. 

Blade  displacement  in  the  axial  direction,  ft. 

Blade  displacement  in  the  lead  lag  direction,  ft. 

Helicopter  forward  speed,  ft/sec. 

Blade  displacement  in  the  flap-wise  direction,  ft. 

Helicopter  gross  weight,  lb. 

Blade  section  angle  of  attack,  rad. 

Longitudinal  shaft  tilt  relative  to  wind  axis,  rad. 

Rotor  coning  angle,  rad. 

Lateral  and  longitudinal  rotor  disk  tilt  angle,  respectively,  rad. 
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Oic^Ou  Lateral  and  longitudinal  cyclic  trim  inputs,  respectively,  rad. 

0  75  Collective  blade  pitch  at  75%  radius,  rad. 

/i  Advance  ratio,  V/QR. 

(T  Rotor  solidity  ratio,  NbCm/T^R- 

^  Blade  twist,  rad/ft. 

<j)^  Lateral  shaft  tilt,  rad. 

p  Air  density,  slug/ ft^. 

Rotor  rotational  speed,  rad/sec. 

Tail  rotor  collective  control  setting,  rad. 


III.  Formulation  and  Analysis 


The  analysis  of  a  bearingless  rotor  system  is  more  involved  than  that  of  a  hingeless  or  an  articulated 
rotor.  A  rigid  rotor  blade  model  with  equivalent  hinge  offset  and  spring,  which  has  been  widely  used  for 
hingeless  and  articulated  helicopter  rotor  system  analyses,  becomes  less  satisfactory,  even  for  approximate 
analysis  of  the  bearingless  rotor.  This  is  because  of  the  multiple  load  paths  near  the  blade  root  and  non¬ 
linear  bending-torsion  couplings  resulting  from  the  large  twisting  motion  of  the  flexbeam.  In  particular, 
the  bending-torsion  couplings  are  very  sensitive  to  pitch-control  configuration  parameters  and  can  also  vary 
significantly  with  changes  in  the  operating  conditions.  Therefore,  all  the  important  structural  elements  of 
the  bearingless  rotor  need  to  be  accurately  modeled  to  capture  the  load  redundancies  and  nonlinear  couplings. 


In  the  present  study,  the  UMARC  comprehensive  aeroelastic  analysis  code  [1]  is  used  to  investigate 
the  effects  of  hub  flexbeam  damage  on  the  dynamic  behavior  of  the  bearingless  main  rotor  system.  The 
finite  element/blade  element  formulation  of  the  UMARC  code  permits  the  analysis  of  bearingless  rotor  with 
complex  hub  configurations  involving  redundant  load  paths.  Recent  applications  of  the  UMARC  code  for 
analyzing  the  dynamics  and  aeroelastic  stability  of  bearingless  rotor  helicopters  and  correlations  with  test 
data  are  available  in  References  3  and  4.  Here,  our  main  efforts  are  1)  to  investigate  the  dynamic  behavior 
of  the  bearingless  rotor  system  with  a  ballistically  damaged  flexbeam  using  the  UMARC  code  as  the  engi¬ 
neering  analysis  tool,  and  2)  to  evaluate  the  results  in  the  context  of  ballistic  vulnerability  assessment. 
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1.  Formulation 


The  formulation  for  the  blade  and  fuselage  equations  of  motion  is  based  on  Hamilton's  principle  [1]. 


/  (6U-6T--6W)dt  =  0  (1) 

Jti 

in  which  SU,  6T  and  6W  are,  respectively,  the  variations  in  the  strain  energy,  the  kinetic  energy  and  the 

virtual  work  done  by  external  forces.  The  expressions  for  6U  and  6T  are  given  in  Reference  1.  The  virtual 
work  6W  can  be  expressed  as 

rR 

8W  —  /  {L^6u^  L^bv-\-  L^6w  +  Mf6<l>)dx  (2) 

Jo 

in  which  and  are  the  external  aerodynamic  loads  distributed  along  the  length  of  the  blade 

in  the  axial,  lead  lag,  flap  and  torsion  directions,  respectively.  In  the  present  study,  aerodynamic  loads  are 
calculated  using  quasi-steady  strip  theory.  A  linear  wake  model  proposed  by  Drees  [5]  is  used  to  calculate 
the  rotor  inflow  distributions.  Dynamic  stall  models  are  also  included  to  capture  some  of  the  unsteady 
aerodynamic  effects. 


2.  Finite  Element  Discretization 

This  section  describes  the  finite  element  method  used  in  the  derivation  of  the  blade-governing  equa¬ 
tions.  In  particular,  the  finite  element  analysis  is  used  here  to  remove  the  spatial  dependence  from  the 
(blade-governing)  equations.  To  achieve  this,  the  bearingless  rotor  blades,  including  flexbeams,  torque  tubes, 
and  main  blades,  are  discretized  into  a  number  of  beam  elements. 


From  the  Hamiltonian  formulation,  the  integrand  of  Equation  (1)  can  be  expressed  as 


A  =  6U-6T-6W 


(3) 


And  in  discretized  form 


A  =  ^A, 

»  =  1 


(4) 


4 


in  which 


A,-  :=  6Ui  -  6Ti  -  SWi  (5) 

and  SUi^STiySWi  are,  respectively,  the  strain  energy,  kinetic  energy,  and  virtual  work  contribution  of  the  ith 
element,  and  n  is  the  number  of  beam  elements. 


Each  beam  element  has  15  degrees  of  freedom  and  consists  of  two  end  nodes  and  three  internal  nodes 
(see  Figure  2).  Each  end  node  has  six  degrees  of  freedom  which  are  u,  t;,  v',  ly,  w\  two  internal  nodes  for 
the  u  (axial)  degree  of  freedom;  the  remaining  internal  node  is  for  the  (f)  (torsion)  degree  of  freedom.  At  the 
element  boundary,  there  is  a  continuity  of  displacement  for  the  axial,  lag,  and  flap  deflections  and  twist,  and 
a  continuity  of  slope  for  the  flap  and  lag  bending. 


In  the  formulation,  the  deflections  at  any  point  within  an  element  can  be  expressed  in  terms  of  the 
shape  functions  and  the  nodal  degrees  of  freedom.  For  the  ith  element 


•  (6) 

. ««)  J 

in  which  the  nodal  vector  qe  consists  of  the  15  degrees  of  freedom  for  the  beam  element,  and 


qe^  =  [‘Ul,U2,U3,U4,Vi,v'i,V2,V2,Wi,w[,W2,w'.2,^l,^2,^3\ 


(7) 


The  shape  function  matrix  H  consists  of  the  shape  functions  for  the  beam  element  degrees  of  freedom 


and  has  the  form 


H  = 


H^is)  0  0  0 

0  H^{s)  0  0 

0  0  Hlis)  0 

0  0  0  Hjis) 


in  which  the  shape  function  vectors  are 


(8) 
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H4s)=< 


> 


(9) 


—4.5s®  +  9s®  —  5.5s  +  1 
13.5s®  -  22.5s®  +  9s 
—13.5s®  +  18s®  —  4.5s 
4.5s®  —  4.5s®  +  s 


H^(s)  =  H^{s)  =  { 


2s®  -  3s®  +  1 
/,(s®  -  2s®  +  s) 
-2s®  +  3s® 
/,(s®-s®) 


{2s®-3s  +  l  I 
-4s®  +  4s  \  . 
2s® -s  J 


The  local  spatial  coordinate  s  of  the  ith  element  is  given  by 


(10) 


(11) 


s  = 


h 


in  which  /,•  is  the  length  of  the  ith  beam  element;  Xi  is  the  local  coordinate  of  the  beam  element,  and 
0  <  X,-  <  The  shape  functions  for  the  lag  and  the  flap-bending  deflection  are  the  Hermitian  polynomials 
which  allow  the  deflections  and  slopes  for  the  corresponding  degrees  of  freedom  to  be  continuous  across  the 
beam  element.  Since  only  the  continuity  of  displacement  is  required  for  axial  deflection  and  twist,  Lagrangian 
polynomials  are  used  as  shape  functions  for  these  degrees  of  freedom. 


To  preserve  the  symmetry  in  the  matrices,  the  virtual  displacements  6u,6vj6wj  and  6^  over  the  ith 
element  are  expressed  in  the  same  set  of  shape  functions  as 


61/(5) 

6i/;(5) 


}  =  H(s)  ^qe 


in  which  the  virtual  nodal  vector  6qe  is  the  variation  of  qe* 


(12) 


Substitution  of  Equations  (6)  and  (12)  into  Equation  (5)  gives 


Ai  =  A,(qe,qe,qe,<5qe)  •  (13) 

Following  Equation  (4),  the  assembly  of  the  beam  elements  yields  a  set  of  nonlinear  ordinary  differential 
equations  in  terms  of  the  nodal  displacements.  Let  q  denotes  the  global  displacements  for  the  rotor  blade 
and  6q  the  associated  global  virtual  displacements.  The  global  equations  can  be  expressed  implicitly  as 
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A  =  A(q,q,q,(5q)  . 


(14) 


Derivation  of  the  blade-governing  equations  is  completed  by  applying  the  geometric  boundary  condi¬ 
tions  to  Equation  (14).  The  rotor  configuration  often  dictates  the  boundary  conditions.  In  case  of  bearingless 
rotors,  three  sets  of  special  boundary  conditions  are  required: 

1)  Constraints  at  the  in-board  end  of  the  torque  tube. 

2)  Compatibility  conditions  at  the  junction  (clevis)  between  the  flexbeam  and  the  main  blade. 

3)  Root  boundary  conditions  (cantilevered)  at  the  in-board  end  of  the  flexbeam. 

The  in-board  end  of  the  torque  tube  is  attached  to  the  pitch  link.  The  in-board  end  of  the  main  blade  is 
attached  to  the  outboard  end  of  the  flexbeam  and  torque  tube  through  a  rigid  clevis.  This  implies  that  there 
is  continuity  of  the  displacement  and  bending  slopes  between  two  sections  of  the  blade.  The  cantilevered 
conditions  at  the  root  end  of  the  flexbeam  imply  that  the  axial  deflection,  flap-bending  deflection  and  slope, 
lag-bending  deflection  and  slope,  as  well  as  elastic  torsion  equal  zero. 


After  applying  the  boundary  conditions,  Equation  (14)  takes  the  form 


A  =  ^q’’(MGq+CGq+KGq-FG)  (15) 

in  which  Mq,Cq,Kq  are  the  global  mass,  damping,  and  stiffness  matrices,  respectively,  and  Fq  is  the 
global  force  vector.  Substituting  the  above  equation  into  Equation  (1)  and  realizing  that  the  blade  steady 
response  is  periodic  for  a  rotor  operating  in  a  steady  level  flight  gives 


/  "  «q(MGq  +  CGq  +  KGq  -  Fg)  =  0  .  (16) 

Jo 

Note  that  the  above  equation  has  been  expressed  in  terms  of  the  rotor  temporal  coordinate  $ . 


3.  Normal  Mode  Equations 

To  reduce  computational  time,  the  blade  finite  element  equations  are  transformed  into  the  normal 
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mode  space  using  a  few  coupled  natural  vibration  characteristics.  The  blade  natural  vibration  modes  are 
obtained  by  solving  an  algebraic  eigenvalue  problem.  By  neglecting  the  external  loads  and  the  damping 
matrix,  the  blade-governing  equations  can  be  expressed  as 


M^q  +  K|.q  =  0  .  (17) 

The  solution  of  (17)  is  expressed  as 


Then  Equation  (17)  can  be  equivalent  to  the  following  algebraic  eigenvalue  problem 

=  (M^)"^Kf5q  (18) 

in  which  the  superscript  s  denotes  structural  quantities.  These  matrices  are  the  inertial  and  stiffness  matrices 
and  are  symmetrical.  Therefore,  this  eigenvalue  problem  can  be  solved  using  Jacobi’s  method  [6].  The 
resulting  eigenvalues  are  real,  and  the  eigenvectors  are  also  real  and  orthogonal.  These  eigenvectors  represent 
the  blade  natural  modes,  and  the  square  root  of  the  associated  eigenvalues  represents  the  blade  natural 
frequencies. 

To  apply  the  normal  mode  transformation,  the  blade  global  displacement  vector  q  is  related  to  the 
modal  displacement  vector  p  by  the  normal  mode  transformation 


q  =  $  p  (19) 

in  which  $  is  an  A^g  x  matrix  formed  with  the  m  blade  free  vibration  modes.  Applying  the  normal  mode 
transformation  to  Equation  (16)  results  in 


in  which 


6p  (Mp 


+  Cp  +  Kp  -  F)  d«'  =  0 


(20) 


M  = 
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c  = 

K  = 

F  =  $^Fg  (21) 

are,  respectively,  the  modal  mass,  damping,  stiffness  matrices  and  load  vector.  In  a  forward  flight  condition, 
M,  C,  and  K  contain  periodically  time-varying  coefficients,  i.e., 


M(^)  =  M(^-h27r) 

C(^)  =  C(^  +  27r) 

K(^)  =  K(^  +  27r)  , 

For  convenience,  all  nonlinear  terms  are  put  into  the  load  vector  F. 


4.  Finite  Element  Discretization  in  Time 


In  forward  flight,  the  blade  modal  equations  are  coupled,  nonlinear  and  contain  periodic  terms.  A  finite 
element  in  time  method  based  on  Hamilton’s  weak  principle  is  formulated  to  calculate  steady  response.  The 
time  period  for  one  revolution  (i.e.,  2  tt  radians)  is  discretized  into  a  number  of  time  elements,  and  within 
each  element  there  is  a  polynomial  distribution  for  response.  From  the  assembly  of  elemental  properties,  the 
periodicity  of  response  is  imposed  by  connecting  the  first  and  last  elements.  This  yields  a  set  of  nonlinear 
algebraic  equations  that  are  solved  numerically  using  a  Newton- Rap hson  procedure.  This  method  is  efficient 
and  robust,  especially  for  nonlinear  periodic  systems. 

The  blade  modal  equations.  Equation  (20),  are  integrated  by  parts  to  yield 

For  a  periodic  solution,  the  right-hand  side  of  the  above  equations  is  zero, 
as 


(22) 


and  they  can  be  expressed 


A2ir 

Jo  ‘ 


=  0 


(23) 
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in  which 


y  = 


and 


Using  the  finite  element  method,  the  time  interval  of  one  rotor  revolution  of  2'ir  is  divided  into  a  number 
of  time  elements.  It  follows  that  the  blade-governing  equations  are  expressed  as  the  sum  of  the  elemental 
expressions  for  all  the  time  elements.  In  particular,  Equation  (23)  takes  the  form 


V  /  6yf  Q<  d’J  =  0  (25) 

in  which  =  0,  ’^jv.+i  =  27r,  and  Ne  is  the  number  of  time  elements  being  used.  The  first  order  Taylor 
series  expansion  of  Equation  (25)  about  a  steady  state  value  of  y©  =  [Po  Po  J 


hj  Qt(yo  +  Ay)  = 


hj  [Qt(yo)  +  Kt,(yo)Ay]  d®  =  0 


(26) 


in  which 


K,,= 


ep 


4Sr-K 


eF  _  p 
ep 
M 


(27) 


For  the  ith  time  element,  the  temporal  variation  of  the  modal  displacement  vector  can  be  expressed  in 
terms  of  the  shape  functions  and  the  temporal  nodal  displacement  vector  as 


in  which  the  local  temporal  coordinate  s  for  the  ith  time  element  is 

^  ~  ’j'.+i  -  ’j'i 


(28) 
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and  0  <  s  <  1.  Note  that  '^,+i  -  $,•  is  the  time  span  of  the  ith  time  element.  In  Equation  (28),  H(s)  is  the 
temporal  shape  function  matrix  which  has  the  form 

H=  [.ffilmi  •  •  • ,  .fffij+llmj 

in  which  Im  is  the  m  x  m  identity  matrix,  and  m  is  the  dimension  of  the  modal  displacement  vector.  In 
the  above  expression,  the  subscript  n*  refers  to  the  order  of  the  polynomial  used  in  approximating  the 
temporal  variation  of  the  modal  displacement  vector.  In  particular,  if  an  nth  order  polynomial  is  used  in  the 
approximation,  then  ti  +  1  nodes  per  degree  of  freedom  are  required  to  completely  describe  the  variation  of 
pj(’®)  within  the  time  element.  Therefore,  H  is  an  m  x  m(nt  +  1)  matrix,  and  ii  is  an  m(m  +  1)  vector. 

In  the  present  analysis,  fifth  order  polynomials  are  used  in  the  approximation,  and  this  results  in  a 
total  of  six  nodes  for  each  time  element  -  two  end  nodes  and  four  internal  nodes.  The  shape  functions  fall 
into  the  Lagrangian  polynomial  group  and  are  expressed  sis 

Hi(s)  =  (-625s®  +  18755“*  -  21255^  +  11255^  -  274s  +  24)  /24 

Nils)  =  (3125s®  -  87505^  +  8875s®  -  3850s®  +  600s)  /24 

Hals)  =  (-3125s®  +  8125s^  -  7375s®  +  2675s®  -  300s)  /12 

HaIs)  =  (3125s®  -  75005“*  +  6125s®  -  1950s®  +  200s)  /12 

/fsjs)  =  (-3125s®  +  6875s^  -  5125s®  +  15255®  -  150s)  /24 
/76(s)=  (625s® -  1250s^  +  875s®  -  250s®  +  24s) /24  . 


With  this  description,  the  modal  velocity  vector  can  be  expressed  as 


P.-  =  H(Vf)^i 


(29) 


In  a  similar  feishion,  the  expressions  for  the  virtual  modal  displacement  and  velocity  vector  are  expressed 


as 


6pi  =  H('^)«e, 

5p.- =  H($)5^..  . 

Substituting  Equations  (28)  to  (31)  into  Equation  (25)  yields 


N’’  [Q.  +  K„-N]  Aiid^ 


0 


(30) 

(31) 


(32) 
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in  which 


For  the  periodic  solution,  the  boundary  condition  for  the  temporal  finite  element  equations  is 


p(0)  =  p(27r)  . 


Finally,  since  the  S^i  are  arbitrary  foi  i  =  1, . . . ,  Ne,  Equation  (32)  takes  the  form 


in  which 


+  Kp  =  0 


■5/.. 


,N 


.=1 


(36) 

1  =  1 

Equation  (36)  is  a  set  of  nonlinear  algebraic  equations  and  can  be  solved  using  the  Newton’s  method. 
In  particular,  are  solved  iteratively,  and  the  subsequent  solutions  are  obtained  using 


During  the  iteration  process,  the  vehicle  trim  controls  are  also  updated  according  to  the  coupled  trim 
procedure  (for  details,  see  Reference  7). 
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5.  Aeroelastic  Analysis 


Aeroelastic  analysis  consists  of  two  phases,  vehicle  trim  and  steady  response,  and  is  calculated  as  one 
coupled  solution  using  a  modified  Newton  method. 

Vehicle  Trim 

Propulsive  trim,  which  simulates  the  free  flight  condition  of  the  vehicle,  is  used  to  calculate  rotor  con¬ 
trol  settings.  The  solution  is  determined  from  vehicle  overall  equilibrium  equations:  three  force  (vertical, 
longitudinal,  and  lateral),  and  three  moment  (pitch,  roll  and  yaw)  equations.  For  a  specified  weight  coeffi¬ 
cient  (Cti;)  and  advance  ratio  (^),  the  trim  solution  calculates  the  shaft  tilt  angles  (a^,!^,),  the  pitch  control 
settings  (^o,  ^i«)}  and  the  tail  rotor  collective  control  setting  {Btr)-  These  trim  values  are  calculated 

iteratively  using  the  modified  hub  forces  and  moments  including  the  blade  elastic  responses. 

Steady  Response 

The  steady  response  involves  the  determination  of  time-dependent  blade  positions  at  different  azimuth 
locations  for  one  rotor  revolution.  To  reduce  computational  time,  the  finite  element  equations  are  trans¬ 
formed  into  normal  mode  equations,  based  on  the  coupled  natural  vibration  characteristics  of  the  blade,  as 
described  in  Section  3.3.  These  nonlinear  periodic  coupled  equations  are  solved  for  steady  response  using  a 
finite  element  in  time  procedure  based  on  Hamilton's  principle  in  weak  form.  One  rotor  revolution  is  divided 
into  a  number  of  azimuthal  elements.  The  assembly  of  elements  results  in  nonlinear  algebraic  equations, 
which  are  solved  using  the  Newton- Raphson  procedure  as  described  in  Section  3.4. 

6.  Analysis  of  Damaged  Bearingless  Rotor  System 

Most  of  the  comprehensive  aeroelastic  codes  developed  to  calculate  helicopter  dynamic  response  and 
stability  employ  the  assumption  that  the  response  of  all  the  blades  is  identical  with  an  appropriate  phase 
shift  for  each  blade.  Thus,  in  the  rotor  response  calculation,  a  set  of  coupled  flap-lag-torsion  equations 
corresponding  to  a  single  blade  is  used.  In  the  study  of  ballistically  damaged  rotor  system,  the  motion  of 
each  blade  has  to  be  represented  by  an  independent  set  of  equations,  and  their  response  calculated  individ¬ 
ually.  Otherwise,  individual  blade  can  not  be  represented.  This  type  of  solution  approach  is  refered  as  a 
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“multi-blade  formulation.”  This  formulation  also  drastically  increases  the  size  of  equations  to  be  solved. 
For  example,  the  size  of  the  system  equations  quadruples  for  a  four-bladed  rotor.  In  the  present  study,  a 
multi-blade  formulation  developed  in  Reference  8  is  used  to  calculate  the  rotor  hub  forces  and  moments. 
Therefore,  the  motions  of  the  individual  blades  are  considered  in  the  calculation  of  hub  loads  and  moments. 


Flexbeam  Ballistic  Damage  Model 


Ballistic  damage  to  the  hub  flexbeam  is  simulated  as  a  reduction  in  structural  capabilities  of  the  af¬ 
fected  blade.  This  is  done  by  reducing  the  mass,  bending,  and  torsional  stiffnesses  (m^,  Ely^Eh^  and 
GJ)  of  damaged  flexbeam  section.  These  values  are  estimated  based  on  the  prescribed  damaged  configu¬ 
ration,  either  observed  from  ballistic  tests  or  by  an  analytical  method,  and  used  as  input  to  the  UMARC 
code.  Since  the  flexbeam  does  not  contribute  to  the  generation  of  aerodynamic  forces,  aerodynamic  effects 
of  damage  are  neglected.  However,  effects  of  flexbeam  damage  on  blade  aerodynamic  loading  distributions 
(i.e.,  aeroelasticity)  are  included  in  the  UMARC  analysis.  The  composite  damage  model  used  in  the  present 
study  is  constructed  in  a  global  sense.  It  is  based  on  the  assumption  that  ballistic  damage  to  the  composite 
flexbeam,  including  laminate  delaminations,  fiber  breakages,  matrix  transverse  crackings,  and  interactional 
fiber-matrix  failures,  can  be  collectively  represented  as  stiffness  reductions  as  well  as  material  discontinuities. 
More  elaborate  treatments  of  individual  failure  mode  of  composite  materials  are  available  in  References  9 
and  10.  Here,  our  main  efforts  are  to  represent  the  level  of  degradation  in  the  composite  flexbeam  component 
attributable  to  ballistic  damage  and  to  investigate  the  effects  on  overall  rotor  and  helicopter’s  performance 
(i.e.,  02,3  mapping). 


7.  Solution  Procedure 

The  following  procedure  was  used  to  study  the  effects  of  flexbeam  ballistic  damage  on  bearingless  rotor 
blade  dynamic  response,  hub  loads,  and  moments  in  forward  flight  conditions. 


1)  With  prescribed  input  data,  vehicle  trim  equations  are  calculated  using  rigid  blade  flapping  (a  starting 
point). 


2) 


Using  control  inputs  from  the  vehicle  trim  solution 


of  Step  1)  and  prescribed  blade  damage  condition, 
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the  blade  modal  information  and  rotor  nonlinear  steady  responses  are  calculated.  The  results  give  detailed 
individual  blade  responses  at  different  span-wise  and  azimuthal  positions. 


3)  Hub  loads  and  moments  are  calculated  using  elastic  rotor  responses.  Then,  the  vehicle  trim  values  and 
blade  responses  are  recalculated  iteratively  using  the  modified  hub  forces  and  moments.  This  step  is  repeated 
until  a  converged  solution  is  obtained. 


4)  Steps  1)  to  3)  are  repeated  for  different  damage  configurations. 


More  details  about  analysis  and  results  of  several  case  studies,  including  main  rotor  blade  damages, 
are  reported  in  References  11  and  12. 


IV.  Results  and  Discussion 


Numerical  results  are  obtained  for  a  five-bladed  soft  in-plane  bearingless  rotor  system.  Results  are  first 
calculated  for  a  undamaged  (baseline)  rotor  configuration.  Results  then  are  calculated  for  this  rotor  with 
simulated  flexbeam  damage  conditions,  and  the  effects  of  damages  are  assessed.  Some  important  parameters 
of  the  rotor  used  in  the  present  study  are  given  in  Table  1. 


For  the  calculation  of  dynamic  response,  each  rotor  blade  is  discretized  into  16  beam  elements:  6 
for  flexbeam,  5  for  torque  tube,  and  5  for  main  rotor  blade  (see  Figure  3).  For  normal  mode  reduction, 
six  coupled  rotating  natural  modes  are  used,  primarily  comprised  of  three  flaps,  two  lags,  and  one  torsion 
mode.  For  periodic  response,  one  cycle  of  time  is  discretized  into  eight  time  elements,  and  each  time  element 
represents  a  fifth  order  polynomial  distribution  of  motion. 
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Table  1:  Bearingless  Rotor  Helicopter  Characteristics. 


Number  of  blades,  Nt 

5 

Root  cutout,  Vc/R 

0.3 

Mean  blade  chord,  c/R 

0.064 

Solidity,  a 

0.102 

Lock  number,  7 

5.09 

Thrust  ratio,  Ct/o' 

0.075 

Rotational  speed,  Q 

355  rpm 

Rotating  flap  natural  frequency 

1.09/rev 

Rotating  lag  natural  frequency 

0.65/rev 

Rotating  torsion  natural  frequency 

5.74/rev 

Table  2:  BMR  Blade  Natural  Frequencies  (per  rev). 


Mode  Number 

Baseline 

Damage  1 

Damage  2 

1  (lag) 

0.652 

2  (flap) 

3  (coupled) 

4  (coupled) 

2.934 

3.051 

5  (coupled) 

5.056 

4.897 

4.299 

6  (torsion) 

5.741 

5.907 

4.998 

In  the  present  study,  two  different  flexbeam  damage  conditions  are  considered.  Damage  1  represents  a 
condition  in  which  the  flexbeam  suffers  25%  reduction  in  structural  capability.  Damage  2  represents  a  more 
severe  damage  condition,  in  which  the  bending  and  torsional  stiffness  of  the  flexbeam  element  is  further 
reduced  to  50%  of  the  pre-damage  values. 


Effects  of  flexbeam  damage  on  the  mode  shapes  of  the  bearingless  main  rotor  blades  are  shown  in 
Figures  4  through  12.  Since  these  modes  are  coupled  ones,  each  flap,  lag,  and  torsional  component  is  shown. 
For  comparison  purposes,  the  results  of  the  baseline  (undamaged)  case  are  also  shown  along  with  the  two 
damage  cases.  It  is  shown  that  the  flexbeam  damage  affects  all  the  modes  of  the  blade.  In  particular,  the 
torsional  mode  is  greatly  affected  by  the  damage  (see  Figures  11  and  12).  This  is  because  of  a  large  coupling 
effect  between  the  bending  and  the  torsional  modes  attributable  to  the  flexbeam.  The  natural  frequencies 
of  the  baseline  and  two  damaged  cases  summarized  in  Table  2. 


From  these  mode  shapes  plots,  one  observes  that  all  components  of  blade  dynamic  motions,  including 
flap,  lag  and  torsion,  are  strongly  coupled.  It  is  also  interesting  to  note  that  this  coupling  phenomenon,  usu¬ 
ally  observed  in  higher  structural  modes,  is,  with  damage,  also  present  in  the  fundamental  mode  shapes.  For 


example,  significant  torsional  components  are  present  in  the  flap  modes.  This  is  thought  to  be  attributable 
to  the  bearingless  rotor  blade  design,  in  which  the  flexbeam  plays  a  strong  role  in  blade  dynamics. 


The  peak-to-peak  amplitudes  of  structural  bendings  are  shown  in  Figures  13  through  15  for  the  undam¬ 
aged  and  both  damaged  cases.  Flap  bending,  lag  bending,  and  torsional  moments  (nondimensionalized  by 
mQ^R^)  acting  on  the  flexbeam  station  {x/R  =  0.1)  are  shown  for  three  forward  flight  conditions  (advance 
ratio  //  =  0.1,  0.2,  0.3).  These  moments  were  calculated  using  a  force  summation  method.  Therefore,  they 
contained  all  the  aeroelastic  information  of  the  rotor  system,  including  rotor  controls  and  wake  distributions. 
The  damage  effects  on  the  flap-bending  moments  are  shown  to  be  large  as  compared  to  the  lag-bending  and 
torsional  moments.  One  also  observes  that  the  peak-to-peak  amplitudes  of  these  moments  are  increased  as 
the  vehicle's  forward  speed  increases. 
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V.  Summary 


Effects  of  flexbeam  ballistic  damage  on  a  soft  in-plane  bearingless  main  rotor  blade  response  and  loads 
were  investigated  using  the  UMARC  comprehensive  aeroelastic  analysis.  The  analysis  was  based  on  finite  el¬ 
ement  theory  in  space  and  time  coordinates.  The  bearingless  main  rotor  (BMR)  system,  including  flexbeams, 
torque  tubes,  and  main  rotor  blades,  was  modeled  as  a  number  of  elastic  beam  finite  elements,  wherein  each 
beam  element  undergoes  flap  bending,  lag  bending,  elastic  twist,  and  axial  deflections.  Aerodynamic  forces 
on  rotor  blades  were  calculated  using  quasi-steady  aerodynamic  theory  with  a  linear  inflow  model.  Flexbeam 
ballistic  damage  was  simulated  as  the  changes  in  the  span-wise  distribution  of  mass,  bending,  and  torsional 
stiffness.  Results  were  calculated  for  the  baseline  undamaged  rotor  system  and  two  hub  flexbeam  damage 
conditions. 

It  was  shown  that  hub  flexbeam  damage  affected  all  the  modes  of  the  blade,  and  these  effects  were 
more  distinct  in  the  higher  modes.  It  was  also  found  that  the  bending-torsion  coupling  effects  were  quite 
significant  in  the  bearingless  rotor  dynamics. 


18 


) 


1. — Schematic  of  Bearingless  Main  Rotor  (BMR)  System. 
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Local  Degrees  of  Freedom 


Figure  3;  Finite  Element  Model  of  Bearingless  Main  Rotor. 
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Deflections  Deflections 
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First  Lag  Mode  (Baseline:  Undamaged) 


Figure  4:  Lag  Mode  of  Baseline  Blade  (Undamaged). 
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Figure  5:  Lag  Mode  of  Damaged  Flexbeam  Blade  (Case  1). 
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Radial  Station,  x/R 


Figure  6:  Lag  Mode  of  Damaged  Flexbeam  Blade  (Case  2). 


Figure  7:  Flap  Mode  of  Baseline  Blade  (Undamaged). 
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First  Flap  Mode  (Damage  Case  1) 


Figure  8:  Flap  Mode  of  Damaged  Flexbeam  Blade  (Case  1). 


Figure  9:  Flap  Mode  of  Damaged  Flexbeam  Blade  (Case  2). 
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Deflections  Deflections 


Torsional  Mode  (Baseline:  Undamaged) 


Figure  10:  Torsion  Mode  of  Baseline  Blade  (Undamaged). 


Torsional  Mode  (Damage  Case  1) 


Figure  11:  Torsion  Mode  of  Damaged  Flexbeam  Blade  (Case  1). 
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VIBRATORY  FLAP  MOMENTS 


0,4 


Torsional  Mode  (Damage  Case  2) 


Figure  12:  Torsion  Mode  of  Damaged  Flexbeam  Blade  (Case  2). 


Figure  13:  Effects  of  Flexbeam  Damage  on  Flap-Bending  Moments  at  .IR. 
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Figure  14:  Effects  of  Flexbeam  Damage  on  Lag-Bending  Moments  at  .IR 
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Figure  15:  Effects  of  Flexbeam  Damage  on  Torsional  Moments  at  .IR 
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